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NATIONAL

FOREWORD

This Indian Standard which is identical with IS0 1088 : 1985 `Liquid flow measurement in open channels - Velocity-area methods - Collection and processing of data for determination of errors in measurement' issued by the International Organization for Standardization (ISO) was adopted by the Bureau of Indian Standards on the recommendation of Fluid Flow Measurement Sectional Commitee (RVD 1) and approval of the River Valley Division Council. Wherever the words `International Spandard' , Standard' appear referring to this standard, they should be read as `Indian

Comma ( , ) has been used as a decimal a point ( . ) as the decimal marker.

marker,

while in Indian

Standard,

the current

practice

is to use

CROSS

REFERENCES the following International Standards are referred to. Read in their respective places the

In this standard following:

International Standard IS0 748 : 1979 Liquid flow measurement in open channels Velocity area methods IS0 772 : 1988 Liquid flow measurement in open channels Vocabulary and symbols

Indian Standard IS 1192 : 1981 Velocity area methods for measurement of flow of water in open channels (jht revision) IS 1191 : 1971 Glossary of terms and symbols used in connection with the measurement of liquid flow with a free surface (first revision) IS 4890 : 1968 Method for measurement of suspended sediment in open channels

Degree

of Correspondence

-

Identical with some elucidations in Indian Standard Based on earlier IS0 772 version of

-

IS0 4363 : 1993 Liquid flow measurement in open channels Methods for measurement of suspended sediment

Not identical

-

The following International Standards have also been referred to for which there are no corresponding Standards but they may be acceptable for use in conjunction with this standard: IS0 4364 : 1977 IS0 5168 : 1978 ISO/TR 7 178 : I983 Liquid flow measurement of fluid
flow -

Indian

in open channels Estimation

-

Bed material

sampling

Measurement

of uncertainty -

of a flow-rate measurement methods Investigation

Liquid flow measurement of total error

in open channels

Velocity-area
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1

Scope

and field of application

3
3.1

General
Principle

Thii Intemetional Standard rpa&ies a standard basis for tha collection and processingof data for the determination of individual components of the total error in the measurement of liquid flow in open channels bv velocity-areamethods. For determining the discharge in open channels by the velocityarea method, components of the flow need to be measured. The total uncertainty in discharge is a combination of the uncertainties in these components. This International Standard specifies a standard basis for collecting and processing the data required to compute the component uncertainties for determining the total uncertainty in discharge, This International Standard may be used when carrying out an investigation of component uncertainties from data taken from a large sample of rivers in a basin or in a country or for international investigations.

The principle of the velocity-arm mathod consistsin determin-

ing from measurementsthe distributionof the flow velocity in the cross-sectionand tha cross-sectional area, and using these observationsfor the computation of the discharge. The measurementsof the flow velocityare made'in a number of verticals. In each vertical the mean velocity is determined from measurements at a selected number of points. The discharge per unit width may be found by multiplying the mean velocity by the depth in the vertical considered.
Each vertical is assumed to be representative of a segment of the cross-sectional srea. The selection of the number and location of the verticals determines the width of these segments. Assuming that the discharge has remained constant during the measurements, summation of the discharge in the various segments gives the total diicharge through the section.

2

References
measurement in open channak -

IS0 748, Liquid flow Velocitv-area methods.

3.2
in open channels -

Occurrence

of error

IS0 772, Liquid flow measurement Vocabulary and symbols.

When measuring width, depth and flow velocity, errors occur. The application of certain computational methods also introduces errors depending on the assumptions made. A distinction shall be mede between random and systematic errors, resulting ,from the instruments used, the measuring procedures and the processing of data. Random errors are also influenced bv the nature of turbulent flow. The magnitude of random errors can be influenced favourably by the proper selection of instruments and methods. Systematic errors may be constant or variable and they cannot be eliminated by repeating the measurements or by increasing the duration of a measurement. There are, in addition, mistakes due to misreading an instrument or to instrument malfunction.

IS0 4363, Liquid flow measurement in open channels Methods for measurement of suspended sediment. IS0 4364, Liquid flow measurement in open channels material sampling. IS0 5168, Measurement of fluid flow certainty of a flowrate measurement. -

-

Bed

Estimation of un-

ISOlTA 7178, Liquid flow measurement in open channels Velocity-area methods - Investigation of total error.

1
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3.3

Sources

of error (see figure 11

3.4 Determination of the error 3.4.1 Erron in width

of the individual

components

Theoretically the discharge

4=

ss

A w(x,y)dxdy

. . . (1)
The measurement of the width between verticals is normally made by measuring distances from a reference point on the bank. When a tape or tag-line is used, or the movement of the wire attached to a trolley is obsenred, the error depends on the distance but is usually negliiible. Where optical means are used, the errors also depend on the distance measuredbut may be greater. Where the distance is measured by electronic means, a constant error and an error depending on the distance measured occur. The errors are mainly instrument errors. 3.4.2 Errors In depth

where
q

is the unobservabletrue discharge; is the cross-sectionalarea; is the velocity field over width, x, and depth, y.

A

Y fx, y)

In practice, the integral is approxirneted by the summation -m bidi . . . (2)

Q=

z

/El
where Q .is the calculated discharge; b,
d,

In IS0 748, clause 10, a number of sources of error in the measurement of depth is mentioned. Some errors depend on the type and use of the instruments r#i&Such errors are not included in this International

is the width of the i* section; Is the depth

of the P" section;
Errors arise due to the interpolation of the depth between the verticals at which depths are measured. These are included in 3.4.3 cj as "error-tvpe Ill". 3.4.3 Errom in dotermlnatlon of the mean velocity

Vi is the mean velocity in the ith vertical; m is the number of sections.

The error in Q is due to a) arrcre in the measurement of the quantities b,, diand of the individual measurements of the flow velocity necessary for the determination of ii), and bj the approximation of the integral (1) by the summation (2).

These errors consist of three components : al The error due to the restricted measuring time of the local point-velocity in each vertical. Because of turbulence the velocity fluctuates continuously over the wet crosssection. The mean velocity at any point, determined from a

b;

Figure 1 -

Definition sketch

2
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during a certain time inten'vd, is an approximation of the true mean velocitv at that particularpoint. In this International Standard, errors of thb nature ars derrd to as "error-type l"V

~uremmt

where

X'sdis the

percentage uncertainty due to random sampling error of the depth profile (error-type Ill);

b) The error arising from the use of a limited number of sampling points in a vertical. Computation of the mean velocity in a vertical as an average or a weighted average of a number of point velocities results in an approximation of the true mean velocity in the vertical considered. In this International Standard, errors of this nature are referred to as "error-type II" 1). 1.9 The error arising from the restricted number of verticals in which velocities are measured. The horizontal velocity profile between two verticals has to be determined bv interpolation which introduces an error. The values of depth d, and the mean velocity 5 in the vertical are used to determine the discharge par unit width and the discharge through the section i. Summation of the discharges through each section according to equation (2) results in an approximation of the true total discharge. Errors of this nature are referred to as "error-type Ill" 1). 3.4.4 Symbols

is the percentagco uncertainty due to the random XShsampling error of the horizontal velocity profile (errortype III); is the percentage uncertainty of the mean velocity due to the random instrumental error; is the p`eresntege uncertainty due to the random fluctuation in velocity error (error-type I); is the percentage uncertaintv due to the random sampling error of the mean velocity in the vertical (error-type II); is the percentage uncertainty due to the random instrumental error determining the depth of section;

xbt
bi di

is the percentage uncertainty due to the random instrumental error determining the width of section; is the width of section i;

The symbols used are given in the table below.

isthe depth

of section i;

I

Percmtqa uncertalntyl~
Qusntltv random
x'b,

I

rYot.matlc~~
x"b

4
Q

is the mean velocity in section i; is the total discharge.

Width

Mean wlocii

I

x'd,
X.3

I
I

x-d

p7

I
I

The total systematic uncertainty is as follows

: . . . (4)

ll Ail tha uncertsinties used in this International Standard are in terms of percentagesexpressedat the 96 % confidence limits. in this lnternetionsi Standard, X with a subscript refers tc percentage uncertainty,X' refers to percentagerandom uncerteintyand X" refers to percentagesystematic uncertainty. 2) The major source of the systemetic uncertainty in the velocity will arise from errors in the calibrationof the current-meters. 3.5

x

I,

Q

=

f

( yb2 +

yd2

+ x",Z
>

112

The total uncertainty in the discharge is given by

xQ-

f

(

112 + x'b2)

XL2

...

Description

of the uncertainty
If it is assumed that the discharge in the respective sections are almost equal, then equation (31 simplifies with certain approximations (see ISO/TR 7178 and IS0 748) to :

For the determination of the influence on the total uncertainty, the individual components are sufficiently described by their relative mean and relative standard deviations. 3.6 Total uncertainty in discharge

It can be shown (see ISO/TR 7178) that the total percentage random uncertainty in discharge, X'g may be found from the following
where

. . . (6)

X'g is the total random uncertainty in discharge at 95 % confidence level; is the percentage random uncertainty X:, restricted number of verticals used; 11 Error-typesI, II and Ill used in this international Standard have no connection with the statisticei type I and type II errors. due to the
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Xi is the percentage measurement; Xi is the percentege meaaurement;

random uncertainty

in the width

random uncertainty

in the depth

In each vertical this procedure should be carried out at 0.2 0,6 - 0,8 end, where possible, 0,9 times the depth meesured from the surface. Where possible the data should ba obtained during the s8rq 2 000 s period. The meesurementa should be repeated for different discharges.

XL is the percentage random uncertaihty restricted time of exposure used;

due to the The data thus obtained should be indicated in the report form given in annex A. In the case of a continuous recorder, the values at intervals of 10 s should be given, indicating the method of determination.

X' is the percentage random uncertainty due to the reltricted number of velocity points taken in the vertical; XL is the percentage random uncertainty meter rating of velocity. in the current

b

Data on the average

velocity"

NOTE - From equation (31 or M, it csn be ssen that in order to reduce ths total random uncertainty the number of verticals would,nesd to b increassd or an improvement made in the msssuretnsnt of the individuel components, or both.

The average velocity in a vertical may be obtained in various ways. The velocity distribution method is, however, taken as a basis for comparison with the results of other methods generally used or special methods adopted owing to special circumstances. The following procedure is required.

3.7 Evaluation of the error in the individual components
The evaluation of the uncertainty in the indiiual components of the total uncertainty can be obtained by a at8tiatical analysis of a large number of obsewations for a perticular component under operating conditions. Incorporating this procedure into the normal routine measurement is not feasible, therefore centralized proceasing of collected data according to atanderdbed programmes, 8s indicated in this International Standard, is recommended with a view to providing a general standard on the uncertainties of the components within the practical range of measurements.

5.1

Location of the vertical

The vertical chosen for this measurement shall normally be datermined from the known velocity distributions in the geuging croaa-section, so 8a to give velocities which are representative of the whole croaa-section. When the velocity distributions in the gauging cross-section are not known, the vertical taken for this measurement shall be that at maximum depth in the cross-section and at places where the depths are 0.6 and 0.3 times the maximum depth respectively, at the side of the greater segment and not too close to the bank.

4

Data on the local point velocity"

5.2

Distribution

of measuring points

To judge the value of a single velocity measurement the following procedure is required. At each point of measurement on a vertical, an uninterrupted observerion of the velocity over a period of 2 OLJO s, or for a period during which the discharge does not change by more than 5 % of the initial value, whichever is the less, shall be made with a current-meter. Every 10 s a reading of the instrument should be taken, thus giving 200 readings altogether. When pulses are emitted by the current-meter, the number of pulses should be recorded every 10 s; or, when the time is measured at a fixed number of pulses, this time interval should average 10 s. When a continuous record is produced, the complete record should be given and the response characteristics of the electronic instrument stated. The verticals to be taken for this measurement should be the vertical situated at the deepest point and the verticals situated at places where the depths are 0,6 and 0.3 times the maximum depth, both located on the side of the greater segment of the width from the deepest point.

The velocity should be measured at the following points in the vertical : 1) 2) 3) 4) immediately below the surface; at 0,2 times the depth; at 0,3 times the depth; at 0,4 times the depth; at 0.5 times the depth; at 0.6 times the depth; at 0.7 times the depth; at 05 times the depth; at 0.9 times the depth;

5)
6) 71

8)
91

10) near the bed.

II

Reference may also be made to IS0 746 and ISO/TR 7178.

4
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In channels containing weed growth great care shall be taken to ensure that measurements made in the vicinity of the bad ere not affected by weed fouling the current-meter.

6.2

Velocity-distribution

method

5.3 Period of measurement of local point velocities
The period of measurement of local point velocity et any point should be 60 s, or the number of pulses should be that observed in 60 s at 0,6 times the depth.

In this method, the normal procedure for discharge measurement may be used provided the velocity-distribution method or integration method is used for the determination of the average velocity in the vertical. Readings shall be taken every 60 s, from a reference currentmeter at a fixed point, preferably at 0,6 times the depth in the vertical at maximum depth. In addition to the date on the depth obtained by the normal discharge measurement, a continuous profile of the crosssection at the measuring site shall be provided, as indicated in 6.1.

.

5.4

Number of measurements

The measurements in each of the verticals should be made at least five times, preferably consecutively. Measurements affected by navigation should be indicated. These sets of observations discharges. should be made for various

6.3

Praaentation of data

5.5

Presentation of data

The mean velocity should be determined with the use of a planimeter from an adequately large graphical plot (preferably not less than 300 cm2). The type and accuracy of the planimeter should be given, together with the scale of the discharge. The accuracy of the graph paper should be checked. The velocity profiles should be drawn to a scale in such a way that the maximum velocity and the depth are represented by 0,lO m and 0,20 m respectively. Compilation of data should be made in the form shown in annex 6.

For the compilation of data, the form shown in annex C should be used. Correction factors in the table on velocity at the reference point can best be based on the everage value of velocity at the reference point. In this table, the factors are set es e function of time. To obtain the corrected velocity in the table "Mean velocity at verticals", the velocity column shall be multiplied by this correction factor. A graphical representation of the cross-section shall be drawn to an adequate scale; the width of the river on the drawing shall be not less than 0,5 m. The represantation shell indicate the numerical values of depth at the measuring points when a rod has been used, and shall show the location of the verticals and of the reference current-meter. A graphical representation of the measured velocity profiles should also be given. This should indicate the numerical values of the velocities at the measuring points. 6.4 General data

6

Date on the velocity-area

method"

There are two possible ways of determining the accuracy of the velocity-area method, one requiring special measurements, the other mainly using routine meesurements. Wherever possible, data for both should be produced.

To facilitate the interpretation of deviations from the normel oattern of the verious errors. relevant information on the geometry and morphology of the river concerned is required, for example a map of scale 1/10 Ooo of the river approximately 50 times the width of the river upstream and downstream of the measuring site.

7
6.1 Measurement at 0.6 times the depth In this method, the continuous profile of the cross-section at the measuring site is required. This can be obtained by echcsounder measurements or by measuring the depth with a rod at intervals in wide rivers of not more than l/50 of the total width. The horizontal velocity distribution shall be observed by taking velocity readings at 0,6 times the depth at intervals in wide rivers of l/50 of the total width. The readings of the currentmeter shall be made over a period of 120 s. In addition, readings shall be taken from a reference currentmeter at a fixed point, preferably at 0,6 times the depth in the vertical at maximum depth. It shall be read every 60 s.

Integration

method

To determine the standard error in the mean velocity in the verticals obtained by the integration method, a sufficient number of measurements (for example, 50) should be carried out at steady stage in three verticals and the results should be tabulated. The verticals to be taken for this measurement should be the vertical situated at the maximum depth and the verticals situated at places where the depths are 0.6 and 0,3 times the maximum depth, both located on the side of the greater segment of the width from the deepest point. The measurements should be repeated for different discharges. Data of a general character can be compiled in a report form similar to that given in annex A.

1)

Reference may also be made to IS0 748.

5
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8

Calibration curves

10.6 Sounding lines (including suspended current-meters) yield errors due to a) penetration into the bed;

In connection with the study of the instrumenterror, calibration curves together with all calibration points should be given, especially data of successive calibratiqns of a representative current-meter with dates and years of calibrationand the intensity of use.

b) deviations from the ideal conditions for which the correction for downstream drift has been calculated;

cl

shape and suspensionpoint of the lead.

9

Distance measurements

10.6 Echo-soundersyield errors due to a) beam width of the transmitted pulse at the bottom;

No generallyapplicable method of determining the accuracy of distance measurements can be given at present. Detailed descriptionof the method of distance measurement should be given, together with the distances involved, and other relevant factors should be given for theoretical examination. Electronic distance measuring devices give an almost absolutely accurate standard of comparison for distance measure merits. Where these instruments are available, independent research programmes, concerning the uncertainty of different methods of distance measurement, may be carried out and the results stated. The conditions under which the study is carried out should be similar to normal operating conditions in the field.

bl penetration of the pulse into the bed, which is a funclion of the frequency of the pulse and of bed consistency.

11

Date processing

11.1 General The method of data processing for the determination of the total random error in the discharge measurement by velocityarea methods is given. Although the availabilityof computers is assumed, it is posa& to perform the comp&ation process with less advanced maans. Some of these alternatives are indicated. When processing the data, steady-flow conditions are assumed, which means that the true mean value of each of the various quantities remains constant with time. Non-steady trends shall be removed from the data before processing (sea 11.2.21. 11.2 Error-type I 11.2.1 Finite measuring time and distribution of results

10

Depth measurements

The accuracy of depth measurementis dependent on the channel conditions and the method of measurement. In the case of lined chennels the bed conditionsare not likelyto influence the accuracy of the measurement. In natural channels, for example rivers, the configuration of the bed varies in longitudinal as well as transversedirections. In relation to the measuring procedure, it is important to know whether the measurement is carried out from a rigid position or from an anchored launch. In the latter case the influence of the irregularityof the bed may result in a greater contributionto the total error of the depth measurement. Owing to the complex nature of the depth measurements, general directivescannot be given. In carrying out a study, the following considerationsmay give guidance. 10.1 In a riverwith a shifting bed, consecutivemeasurements at one point should be avoided. 10.2 It is advisable to study the bed configuration in the vicinity of the actual measuring point by determining longitudinal and transversesections. 10.3 For all instruments the accuracy of the reading in .relation to the scale intervals should be determined. 10.4 Sounding rods yield errors due to a) b) c)
penetration into the bed; deviation from the vertical position; built-up head due to velocity.

The standard deviation of the fluctuation error due to a finite measuring time is calculated. It is assumed that the means found from the actual measurztments are equal to the hypothetical means over infinite measuring time and that the distribution of the results is of normal (Gaussian) nature. 11.2.2 Correction for non-steady conditions

The mean velocity is calculated from

VAwhere

n

c
i=

n

v,

(7)

1

1;; is the mean velocity during the time interval of measurements; V, n is the observed instantaneous velocity; is the number of observations.

When the velocity V, is plotted against time I,, it can be seen from the graph whether the magnitude of Vi shows a certain

8
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trend which indicates that the conditions during the measure ments were not steady. If so, the observed velocities shall be corrected using vi = rr tti - T, + b

11.2.4

Autocorrelatien

function can be determined from the

The autocorrelation equation

function

9 . * (8)

n-k c B(k) =-& j=' n rvi 1 VP (vi-nwi+,-v,

where ri is the time when the velocity is Vi and in which the constants II and bare determined from the following equations :

. . . (141

. . . (9)
where and p(k)

c i=

is the autocorrelation function; is the time displacement of function;

a=

c vjtr,- c c T)*
7,
i=l

n

n

Viti -

nT~
R

k

n

=

i-1
n

(ti -

c

$2 - l/n
i

c
i= 1

ti

i=l

=l

).
. . . 111)

2

. . . (10)

c
i=l

n

( Vi - VI *

is already known from previous calculations

[see equation (1211.

The observed instantaneous velocity

Vi is thus corrected :

The autocorrelation function is used for the calculation of the standard deviation, as described in 11.2.8. 11.2.6 Effect deviation (1)

Vcorri = Vi -

(Vi -

~7)

where V-, is the corrected velocity assuming a linear trend in the observed velocities tin the following equations Vi denotes Vcon,if correction has been applied).

of measuring

time on standard

11.2.3

8tandard

deviationtr

of velocity

fluctuations

The standard deviation of the velocity fluctuations is calculated by using In

If no computer is available, the influence of the measuring time for the interval kt, in which t, is the initial measuring time and k is an integer, the standard deviation can be determined. For that purpose the mean velocity over an interval kt, is calculated from :
`2i+k1

=

Vi + Vi + 1 + . . . . . + Vi + k _ 1
k
`2i+k-1

2

. . . (15)

SF=

v

E
i=l

(vin-1

vr*
. . . (12)
The magnitudes of

2

are used to calculate the stan-

dard deviation [see equation (1211. procedure can be simplified by using the This method is less common than the procedure described in 11.2.6. 11.2.6 Effect deviation (2) of measuring tlme on standard

The computation equation

n &2_ c i=

i

i='

c
n

"

2
vi

i . . (13)

SF =

1

/
2

n-l

Using the standard deviation S ito) for an initial measuring time t, be 11.2.3) and the autocorrelation function d(k) (see 11.2.4), the standard deviation of the velocity fluctuations can be calculated from :

(= rt ?*t is already known from previous n calculations [see equation (711.

In this equation, the influence of the measuring time on the accuracy of the point velocity is characterized.

1) For the relationshipbetween standard deviation and uncertainty, see IS0 5166. In general, uncertaintiesare expressed at the 95 96 confidence kw8l and usually as percentages. For large values of n, standard deviationsshould be multiplied by 2 and by 100 to obtain percentage uncertainties. The factor is greater than 2 for smell samples.
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11.2.7

Compilation

of results

error for each of the computation formulae are (see also IS0 748) : a) b) c) d) V = vo,a V = 0,5 (LQ + vo,a)

methods.

Some common

The results of (see the table fluence of the velocities by a Uncertainties IS0 5188.

the calculations should be compiled in a table below) which gives an impression of the inmeasuring time on the accuracy of the point comparison of the various standard deviations. should be calculated in accordance with

V = 0,25 vo,2 + 0.5 v,-,,,j + 0,25 v~,~ 1 v = 3 (Vo.2
+ vO,$ + t'o,~,)

The percentage standard deviations are found by dividing the standard deviation by the average value of the variable and multiplying by 2 (95 % level) and by 100 (see IS0 5188).

d

v =

o,l

vsurf

+

0,s

vg.2

+

0.3

v0,6

+ +

0.2

vO 8 + 0.1

v&d

11.3

Error-type

II

f)

V = 0,l v.,rf + 0.2 VO,J + 0.2 v()4 + 0,2 v,,s + 0.1 v,,&

0,2 vo,s +

11.3.1

Approximation

of mean velocity

in the vertical

The relative standard deviation due to the approximation of the mean velocity in the vertical caused by applying a finite number of point velocities is calculated.

The calculated mean velocity is standardized by dividing it by the standard mean velocity obtained according to 11.3.2.

11.3.4 Sampling error due to velocity computation rule

fluctuations

and

11.3.2 Determination the vertical

of the standard

mean velocity

in

As described in clause 5, the standard mean velocity in the vertical is found by determining the area of the vertical velocity profile with the use of a planimeter and taking the depth into account. As a rule the mean velocity can be calculated with the use of the trapezium rule applied on ten point velocities :

From the results of ten measurements under identical conditions of the standardized velocities, the mean and the standard deviation are calculated assuming

. . . (18) i=
1

. . (19) vi = ( vsufi + 2
vo.2 + 2 vo.3 + 2 vo.4 + 2 vo.5 + 2 VO.6 +

i=
where

1

+2VoJ+2Vo,*+2V~,9+

VW)+-

. ..(17)

J Fr,l v,,li

is the number of observations in a series; is the mean relative velocity in a series; is the mean relative velocity in a vertical; standard deviation of the mean relative

where i is a sequence in the series of measurements.

11.3.3

Computation

methods S rBI is the velocities.

By the application of computation rules the mean velocity in the vertical can be approximated by calculation. Comparison of this approximate mean velocity with the standard mean velocity makes it possible to determine the sampling Table Location of measuring point h,,
0.2 63 0.4 0,s 0.6 0.7 0.8 0.9

The sampling error, Sit due to the computation rule is equal to if a series of I sets of measurements are con(1 - P,,,)I sidered.

-

Compilation 0 (2)

of results

XS %

0 (1)

.. ..

6 (k)

x, Ir,)

xs (2 I,J

%

%

x, lkl,) %

8
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The systematic part of the sampling error can be estimated by

For the purpose of determination of error-type discharge is assumed to be the true discharge.

III,

this

f20)

11.4.3

Omieelon

of vertloale

where i, represents the standard deviation of the sampling error due to the velocity fluctuations under steady conditions. The mean standard deviation of all sets of measurements due to velocity fluctuations can be determined from

Omitting verticals from use for the calculation of the discharge, the absolute error can be determined by comparison with the true discharge. The calculation of the discharge shall be based either on the mid-section method or on the mean-section method, as described in IS0 748. As is proved in ISO/TR 7178, the results obtained by usa of these methods show differences which are negligible. The mid-section method, however, shows a time saving and cost saving over tha mean-section method. The nor is standardized by dividing the absolute error by the true discharge.

(211

11.3.6

Sempling

error due to oomputetion

rule

Of a seriesof I sets of measurementsthe standard deviation of the sampling error due to the computation rule can be calculated :

11.4.4

Maan and standard doviationa of error

ss:=&i
i=

(pi1

i,J* -

J

SF2

The r&tive mean error and its standard deviation due to the omieekm of verticals am determined from a number of series of Wts in the same croes-section.

...

(221

(23)
where is the standard deviation of the sampling error due to ttG computation rule; ii is the sampling error in mean velocity in vertical due to the computation rule per series; fi, is the mean sampling error for all skies together; where

S

and

SSkd(m) =

J
i=l

~[y-1]*
I1

'
* . . (24)

S,

is the mean standard deviation of all series together due to velocity fluctuations; is the number of sets of measurements in a series; is the number of measurements in a series.

fi (ml is the relative mean error in the discharge due to the restricted number of verticals; m is the number of verticals considered for the determination of discharge; Qi (ml is the calculated discharge applying m verticals using the ith maasuramant; Qi is the approximated measurement; Sshd(m) true discharge using the ith

I J

11.4

Error-type

III

11.4.1

Restricted

number

of verticals

The standard deviation of the error of the discharge due to measurement in a restricted number of verticals is calculated.

is the relative standard deviation of error due to

error-type Ill, applying m verticals Ssh,r = Ssb* + Ss,`; 11.4.2 True discharge

I

is the number of measurements.

Using velocities, from which non-steady trends have been removed, and the corresponding depths and applying the midsection method, the discharge in a cross-section is calculated according to IS0 748.

The mean error or its standard deviation is plotted against the number of verticals or the maan of groups of verticals considered. From these graphs the relation between the uncertainty and the number of verticals applied may be deduced.

9
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11.45

Critorla applied in choosing the verticals

12
a b

l&t of symbols
coefficient of linear regression; coefficient of linear regression; unobservable true width of section i; unobservable true depth in vertical of section i; number and Ill); of series of measurements (error-types II

The distances between adjacent verticals on a cross-section may be al b) equal; varied, to make the area of all segments equal;

bi
di I

c) varied, to make the discharge through all segments similar; d) varied, to minimize the influence on the discharge of changes in velocity and depth. Items a) and b) need no explanation.

J k

number of measurements per series (error-type II); time displacement interval, etc. ); autocorrelation function (of time

Item cl can be explained as follows: The true total discharge being 100 %, the discharge per segment will be + % when m

number of verticels or sections per cross-section; number of time intervals of measured velocities ferrortype II;

m segments are considered. The cumulative discharge is plotted ageinst the percentage of total discharge, as shown in figure 2, to establish the position of successive verticals. As soon as the location of the ith vertical is fixed, and the number of verticals to be applied and therefore also the percentage of the discharge which should flow through each segment are determined, the location of the (i + 11th vertical may be found. In figure 2 an example is given for a discharge of 8,s. 96. Item d) can be explained as follows: Based on the equation for the total discharge according to the mid-section method

n

Q Qi

total dischsrge; discharge of section i; stochastic sampling error of mean velocity in vertical (error-type III; time i; initial measuring time (basic time interval); velocity at time i or in section i; actual velocity at time i or in section i;

si
ti

to
"i

+
where

brA:AF >

...

c25)

vi

Vcorri actual velocity from which trend is removed; Adi
is the difference in depth in the verticals i and i + 1 respectively; ai width of the ie section; depth of the i* section; calculated mean velocity; unobservable true mean velocity; unobservable velocity; random instrumental error of mean

di
AFi is the difference in velocity in the verticals i and i + 1 respectively. The optimal selection of the verticals as affected by differences in velocity is based on the product bi dj *ii/, and as affected by differences in depth on the product bi vi Adi In the first case the depth in all verticals is used in the calculation of the discharge. In that way the error due to interpolation of the velocity between the verticals can be estimated. In the second case, the error due to interpolation of the depth between the verticals can be determined. The magnitude of the product indicates the influence on the total discharge of differences in velocity and depth respectively. Verticals are omitted, starting with the smallest, depending on the magnitude of the product. V;

vi
6Fi

Fi

total unobservable random fluctuation velocity (error-type I);

error of mean

sv,
L

unobservable random sampling error of mean velocity (error-type II); mean sampling error (error-type II);

&

i (rn) mean error when m verticals are applied (error-D/pe Ill);

10
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p (k) X

autocorrelation function with time diilacernent percentage uncertainty tgeneral); percentage random uncertainty; percentage systematic uncertainty;

k;

x%
'

percentage uncartaintv due to the random sampling error of the horizontal velocity profile (error-type Ill);

, uncwtainw due to the random sampling XSpI lzercemtaga
error of the nwnn velocii in the vertical (error-typell).
NOTE - Duo to the statistical nature of this Intemetional Standard, it is necenary to hava symbols representing observed values and true values of variabk~ The ymbols may therefont not conform to IS0 772. Observationsor resultsof calcubtions using observationsant indicated by a capital letter. Statisti& quantities obtained from obssrvationsare indiited by a small tier with can%. A mean value is indicated bv a bsr.

X' X"
Xh

percentage uncertainty due to the random fluctuation in velocity error ferror-type 1); percentage uncertainty due to the random sampling error of the depth profile &or-type III);

xi&I

11
c
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I
r-This

I

I

segment carries 8.5 % total discharge

10

15
Figure 2a -

20
vd curve

25

30

35

Number of the vertical from initial point

25
Number of the vertical from initial point Figure 2 b -

35

Cumulative Computation

discharge

curve

Figure 2 -

of discharge
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Annex A Local point velocity measuiements (See clause 4)
NOTE -

-

Report form

All values should be converted if necessely to metres and seconds, where appropriate.

Measurements carried out under the responsibility of Address:

. ..,.,.... ..............................

......................................................... ..................................

.......... ..........

Please refer inquiries to Riverconcerned:

:

........,.......................... ... ...

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ..at(l~~onl.....................................

Country:..........................................................

1

General data on the river during
Discharge Area of cross-sectionof stream Average velocity Water surface width Maximum depth Average depth (area/width1

measuring

period ,..... ...... ...... ...... ,..... ,..... .......... .......... .......... .......... .......... .......... .......... .......... . . . m3/s 9. d . . m/s
... ... ... m m m m3/s.m

Slope of stage-dischargecurve (sensitivityin steedy-stage condition1 Water temperature

. .

OC

2

General data on measuring

verticals

A Meximum depth B 0,6 times msximum depth C 0.3 times maximum depth

I

I

I

3

Method
a) b)

Of Velocity

measurement

(P&se

tick as appropriate)

q
cl

Number of pulses, counted every 10 s. Time (average 10 sl for a constant number of pulses, being for A: ._................................. ................................... . . . . . . . . . . . . . . . . . . . . ..,............ pulses. pulses. pulses.

q
c)

B: c:

Continuous velocity record

: : ..... . .. .... . . .

Speed of paper at recorder

.... Response characteristics of the electronic instrumentation :

.mm/s
.. .... . .. . . ...... ...

cl

d)

Integration method.

13
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4

Measuring

apparatus

Type of current-meter (Please give details, such as diameter of propeller and pitch or size of cup, cup distance, and number of cups) :

..................................................................
Serial number : . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
TypeofBuspension:................................................ Date of calibration (Please add rating-curve and equations for calculation of velocity and indicate whether group calibration or direct calibration is intended. Check whether the water temperature during calibration has been indicatedontherating-curve.)......................................................................................... Methodoftimingandassociateduncertainty ............................................................................

5 5.1
5.1.1

Measurements

at verticals

Vertical at maximum depth
General data
Units

0.2
Date of measurement Hour of measurement bcal time) Water-level at start bcal units) Water-level at end flocal units) Water-depth during measurement Depth of instrument under surface Average velocity wer a period of 2ooo s Minimum velocity Maximum velocity Method of velocity measurement : a) or 01 (me clause3 of this annex) Tims interval of readings Method b) Number of pulses Average time interval of readings

0.6

0,s

g,g m m m/s m/s m/s

10

10

10

10

s S

6.1.2

Ducription

of bed conditions
................................................................................... : .................................................................................... load.

Bedmaterial(size,shape,density)l): Sedform(smooth,ripples,dunes) Sediment transport

: yeslno; if yes, type of transport : bed load/suspended

Bed roughness (preferably expressed in Chewy coefficient Cl

: .............................................................

where C V Rh S is the Chezy coefficient; is the mean velocity; is the hydraulic mean depth (area of cross-section of stream over wetted perimeter); is the energy slope.

R~marks:.........................................~........................,................................... ................................................................................................................... 1) See IS0 4363 and 4364.
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6.1.3

Velocity observations
Velocity at 0,2 depth every 10 8. Reading Velocity Serlal No. Reading Velocity Serial No. 101 Reeding Serial No. 151 Reading

_
a)

serlrl
No.

I ' I
'

L I
50

f

L Average v&city

+

Maximum velocity Minimum vel0cit-v

I

Average velocity
Maximum velocity Minimum velccitv I

t t

150

200 bmrege velocity Aaximum velocity Jlinimumvelocity

Average vebcii Maximum v&city I Minimum vebcitv

NOTE - PIeare indicate special circunwtances influencing tha measurements, such as passing of ships.

b)

For measurements at 0,6 -

0.8 and 0.9 depths, give details in the form shown at al above.

5.2

Vertical at 0.5 times the maximum

depth

r
Vslocity

Give details in the form shown under 6.1.

6.3

Vertical at 0.3 times the maximum depth

Givedetaik in the form shown under 5.1.

IS 14573 : 1998 " IS0 1088 : 1985

Annex B Average velocity
NOTE -

measurements (See 5.51

-

Report form

All valuea should he converted if naceaaary to metres and seconds, where appropriate.

Measurements carried out under the responsibility of ........................................... Addrsss: ......................................................... ......................................................... Pleass refer inquiries to Riverconcsrnsd:

: .............................................
.,at(location): ....................................

.................................................

Country: ..........................................................

1

General

data on the river during
Discharge Area of croee-sectionof stream Average velocity Water surface width Maximum depth Average depth fame/width)

measuring

period ............ ............ ............ ............
maIs m2 m/s m m m

Slope of stage-discharge curve (eanaitivity in steady-stage condition) Depth at vertical Average velocity on vertical Distance of vertical from right hank Water temperature

............ ............ ............ ............ ............

ms/s.m m m/s m OC

2

Method
a) b) cl

of velocity

measurement

(Please tick as appropriate)

0
cl 0

Number of pulses counted every SO s. Time (average SO s at 0.6 depth) for a constant number of pulses, i.e. Continuous velocity record : Speed of paper at recorder

...

..

.

pulses.

: .. .. ..... ........ ..

. .. .. : ...

mm/s

Response characteristics of the electronic instrumentation

.

.

.. . .

.

..

cl

d)

Integration method.

18
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3

Measuring

apparatus

Tm of current-meter (Pleaw give detaL, wch (II diimeter of propeller and pitch or size of cup, cup diitance and number of cupe) : . . . . . . . . . . ...* ,.........*
Serial number

........ ... ........ ........ ..... ....... ..

: .................... ........ ....... .. ...... ....... ..

Typeofsuspension:................................................ Date of calibration (Please add rating-curve and equations for calculation of velocity and indicate whether group calibration or direct calibration is intended. Check whether the water temperature during calibration has been indicatedontherating-cuwe.)..................................................................,......................

4 4.1
4.1.1

Measurements

at verticals

Vertical at maximum depth
Generel date 1 Date of measurement Hour of measurement (local time) Water-level at start (local units) Water level at end llocal units) Water-depth during measurement Duration of the local point velocity measurements using method a) or cJ bee cbuw 2 of thb annex1 Number of pulses using method b) based on a duration of 60 s at 0.6 depth of the local point velocity measurements 3 4 6 Unite

80

T t-Y
m s 60 8 60

4.1.2

Dercription

of

bed condition8
...................................................................................

Sedmaterial(size,shape,density)l): Bed form (smooth, ripples, dunes) Sediment transport

: ....................................................................................

: yes/no; if yes, type of transport : bed load/suspended loed. : .............................................................

Bed roughness (preferably expressed in Chezy coefficient, Cl

where C ii Rh S is the Chezy coefficient; is the mean velocity; is the hydraulic mean depth (area of cross-section of stream over wetted perimeter); is the energy slope.

Remarks: .......................................................................................................... ...................................................................................................................

1,

See IS0 4363 and 4364
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4.1.3

Velocity

observations

1 Velocities Reading
1I 2) 3) 4) 5) 6) 7) 8) 9) 10) near surface 0,2 depth 0.3 depth 0.4 depth 0,5 depth 0.6 depth 0.7 depth 0.8 depth 0,s depth near bed21

2 Velocity
Reading Velocity Reading

3
Velocity Reading

4
Velocity Reading

5
Velocity

1f

Average velocity determined graphically 1) 2) I/ \I I/ \I I/ \I I/ \I I/ \I

The observation near the surface has been taken. The observation near the bed has been taken

m below the surface. m above the bed. such as passing of ships.

NOTE -

Please indicate special circumstances influencing the measurements,

4.2

Vertical at 0,8 times the maximum depth

Give details in the form shown under 4.1.

4.3

Vertical at 0,3 times the maximum depth

Give details in the form shown under 4.1.

18
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Annex C Velocity-area method (SW 6.3.) Report form

NOTE - All values should be converted if necessary to metres and seconds, where appropriate.

Measurements carried out undertheresponsibityof ............................................ Address

: ........ ... .. ........ .. ...... ...... ...... ...... ....... ... .........................................................

Please refer inquiries Riverconcerned:

to

: .............................................
.,et(locstion) : ....................................

.................................................

Country: ..........................................................

1

General data on the river during

measuring

period .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... ..........
IAl

Discharge Area of cross-sectionof stream Average velocity Water surface width Maximum depth Average depth (area/width) Slopa of stagedischarge curve (sensitivityin steady-stage condition) Water temperature

N?
m/s m m m

2
0 cl Cl

Method
a) b) c)

of velocity

measurement

(Please tick as appropriate)

Number of pulses counted every 60 or 120 s. Time(average60orl20sat0,6depth)foraconstantnumberofpulses,i.e. Continuous velocity record . . . . . . . . . . . . . . . or .,............. pulses.

: : .. . ., ... . . . . . . . . . . . . . . . , _ . mm/s : . .. . .. ...... .... . ... ......... .. ......

Speed of paper at recorder

Response characteristics of the electronic instrumentation

0

d)

Integration method.

3

Measuring

apparatus

Type of current-meter (Please give details, such as diameter of propeller and pitch or size of cup, cup distance and number of cups) :

.....................
Special number

..I............................__............

: . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ..____...

Typeofsuspension:............................................. Date of calibration (Please add rating-curve and equations for calculation of velocity and indicate whether group calibration or direct calibration is intended. Check whether the water temperature during calibration has been indicatadontherating-curve.).........................................................................................
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4

General data on velocity measurement

4.1 Velocity observations
Unita Date of measurement Water-level at start Meal units) Water-lwel at end kal units) Number of verticals Number of points per verkal Duration of the local point vekcity messurement using method aI or cl (sea clause 2 of this annex) : at 0.6 depth measurement at reference point bv v&&v diition method Number of pulses using method bl, bawd on : duration of 120 s for ths 0,6 dspth messurement duration of 60 s for ths refsmncs point velocity distribution method Profile of cross-ssction messured bv echwowder/rod st intwvb of -

120 60

8 8 S

m

4.2 Description of bed conditions
Bedrnatar&l(elze,&ape,dansity)l~
Bed form (smooth, ripples, dunes) Sediment transport
:`. ..................................................................................

: ....................................................................................
type of transport

: w&no;

if $3,

: bed load/suspended load. : .............................................................

Bed roughness (preferably expressed in Chewy coefficient, C)

c=v
463
where
C V is the Chezy coefficient; is the mean velocity; is the hydraulic mean depth (area of cross-section of stream over wetted perimeter); is the energy slope.

Rh
S

Remerks:.......................................................................................................... ,.....,...,...,.,,.,.,.......,..................................................,..,.,.............................

b

Method of measurement of mean velocity at verticals
0.6 depth measurement. (Duration of individual velocity measurements 120 8.1

(Please tick as appropriate)

cl
0

Velocity-distribution

method.

(Normal practice; please add procedure. 1

1)

See IS0 4Xi

and IS0 4364.
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Mean v&city

at verticals

Serisl No.
1

Time

vuloclty

Corrected
V&City

Depth `1

Serial No.
26

nmo

Valocity

Corrected
VdOCity

Depth `1

25 1) To be filled in when the velocity-distribution method has been applied.

50

Velocity

et reference

point

Timal)

Velocity

Correction facto9

Time

\ieloolty

Correction factor2)

Time

Velocity

Correction factorz)

1) Intervals of e-u5.
2)

. __

Cormctionhas bean made with respect to the velocity at

.

. ...

h.
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